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Abstract 



I present a first-principles theory of diffusion-limited aggregation in two di- 
r^ ' mensions. A renormalized mean-field approximation gives the form of the unsta- 

'^ ■ ble manifold for branch competition, following the method of Halsey and Leibig 

^ ' [Phys. Rev. A 46, 7793 (1992)]. This leads to a result for the cluster dimension- 

ality, D ^ 1.66, which is close to numerically obtained values. In addition, the 

X' 

Vh ' multifractal exponent r(3) = D in this theory, in agreement with a proposed "elec- 

trostatic" scaling law. 



PACS Numbers: 64.60A, 68.70, 05.20 



Diffusion-limited aggregation (DLA) is a model of pattern formation in which 
clusters grow by the accretion of successive random walkers.^ Each random walker 
arrives from infinity, and sticks to the growing cluster at whichever surface point it 
first contacts. Only after the accretion of a walker does the next walker commence 
its approach to the cluster. The clusters thereby obtained are fractal in all dimen- 
sionalities d > 1, and are qualitatively and/or quantitatively similar to patterns 
observed in such diverse phenomena as colloidal aggregation, electrodeposition, vis- 
cous fingering, and dielectric breakdown.^ 

At the heart of the problem of diffusion-limited aggregation is the following 
question: what is the relationship between the scale-invariance of the diffusive 
growth process and the hierarchical structure of the clusters generated by this 
process?"^ A preliminary, and incomplete, answer to this question was provided by 
this author in collaboration with M. Leibig.^ In this work, it was hypothesized that 
the quantitative process by which one branch screens, i.e., takes growth probability 
from, a neighboring branch, has a specific form, independent of the length scale 
on which this process takes place. This assumption allows the development of a 
qualitatively correct theory, which yields multifractal scaling of growth probability, 
as well as agreement with a phenomonological scaling law, the "Turkevich-Scher" 
law, relating the scaling of the maximum growth probability over all sites on the 
cluster to the dimension of the cluster as a whole. ^ 

In this letter, I shall present a more complete and a priori theory of diffusion- 
limited aggregation in two dimensions based upon a specific mean-field calculation 
of the dynamics of branch competition. Because the mean-field approximation is 
implemented on all length scales, it is perhaps better to regard this theory as an 
ansatz solution in the case where certain types of fluctuations on all length scales 
are neglected, while others are included. This specific model allows verification of 



all qualitative aspects of branch competition that were advanced as (reasonable) 
hypotheses in Ref. 4. The result obtained for the dimensionality of the cluster, 
D = 1.66, is within 3% of the oft-quoted value D = 1.71 obtained from the scaling 
of the cluster radius-of-gyration in numerical studies. An additional scaling law (the 
"electrostatic scaling law"), relating the multifractal exponent r(3) of the growth 
measure to the dimensionality D by D = r(3), is seen to be exact within this 
theory.^ 

In the growth process, each particle attaches itself to a unique "parent" particle 
in the pre-existing cluster. Furthermore, the cluster is observed to be a branched 
structure, with no loops and with each particle having asymptotically zero, one or 
two "children", i.e. particles to whom it stands as a parent.^ Very rarely particles 
have more than two children; primarily for reasons of convenience I neglect this 
possibility. 

Consider a particle with two children. Each of the two children separately, 
with all particles descended from each, I term a "branch". Thus these two-child 
particles are parents of two branches, which occupy neighboring regions of space. 
The total number of particles in one branch I term ni, and the total in the other 
77-2. The total number of descendants of the parent particle is thus Ub = ni + n2- 
Now consider the next particle to accrete to the cluster. I say that this particle has 
a total probability pi to stick anywhere on the first branch, and a total probability 
P2 to stick anywhere on the second branch, yielding a total probability pb = pi +P2- 

Let us now consider the normalized quantities x = pi/pb and y = ni/ub- 
Clearly drti/dn = pi, where n is the total number of particles in the cluster, and 
we are neglecting fluctuations of 0(v^^). Thus y obeys the following equation of 
motion: 



x-y. (1) 



dlnub 

The right-hand side of this equation is a function only of x and y. Now x will obey 
an equation of the form 



G{x,y;n; {<!>.}), (2) 



dlnub 

where {(pi} is some parameterization of all of the variables describing the structure 
of the cluster. In ref. 4 we assumed that by averaging the right-hand side of this 
equation over these parameters {(/){}, one obtains dx/dlnub = g{x,y), where the 
right-hand side is now only a function of x and y. Given this function g{x, y), one 
has a closed system of equations describing the evolution of x and y as functions of 
In rife. 

By symmetry, g{x,y) = —g{l — x, 1 — y), so {x,y) = (1/2,1/2) must be a 
fixed point of this process of competition between the two branches. In ref. 4, we 
explored the consequences of assuming that this fixed point is hyperbolic, with the 
unstable manifold emerging from the fixed point terminating in two stable fixed 
points at {x,y) = (0,0) or (x,y) = (1, 1), these latter representing the situation in 
which one branch has been completely screened by the other. This assumption will 
be explicitly verified in the calculation below. 

If the central fixed point at {x,y) = (1/2, 1/2) is hyperbolic, then branch pairs 
which commence their existence (with n^ ~ 1) near the unstable fixed point will 
be quickly drawn onto the unstable manifold. Linearizing the system of equations 
for d{x,y)/dlnnb about the central fixed point, the hyperbolic assumption implies 
that there will be a stable and an unstable direction; the eigenvalue corresponding 
to the latter direction we define to be u. 



When a pair of branches is first created by a tip-sphtting event, its initial 
growth up to the stage at which n?, ^ 1 is determined by comphcated microscopic 
dynamics, which do not recognize the existence of the unstable fixed point. Thus 
we expect the probability that a newly created branch pair will be a distance e^ 
from the unstable fixed point will be p(e ^ l)de oc e^~^de] we are assuming a 
constant probability density of branch creation near the unstable fixed point. This 
assumption has been specifically verified by numerical study in ref. 4. The choice 
of e^ for this initial distance insures that position along the unstable manifold in 
the X — y plane can be parameterized by the variable en^. 

It is possible to relate the eigenvalue u to the cluster dimensionality D by the 
following argument.^ Consider the strongest branch in the cluster, that obtained by 
always following the stronger child (with the larger values of x,y) at each branching. 
The total number of side-branches (or branch points) from such a branch is ~ r, 
where r is the cluster radius. In order that the cluster have a dimension D > 1, a, 
number ~ 1 of these side branches must have a total number of particles ~ n, the 
total number in the cluster. A side branch obeying this criterion must have en ~ 1, 
so that at that branching, both descendant branches are roughly equal in size. The 
probability of this happening at any particular branching is / de p(e) oc n'^, and 
there are ~ r different sidebranchings at which this might occur. Thus rn'^ ~ 1, or 
D = Ijv. 

In order to determine g{x,y), we turn to an explicit description of the growth 
process. ^'^ Suppose that we parameterize the accessible surface of the cluster by 
arc-length s. If a particle attaches at the surface point s', it thereby reduces the 
growth probability at all points s for which |s — s'| > a, where a is the particle size. 
This is because a certain number of the random walks that would have reached s 
previously are now obstructed by the new particle at s' . If the probability that a 
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particle lands at s' is p{s'), and the probability that a random walker goes from s' 
to s without contacting the surface is H{s, s'), this implies that 



dn 



ds' {H{s, s') - h{s)5{s - s')) p^{s'), (3) 



where we have modelled effects on the scale |s — s'| < a by the 5- function, the 
coefficient of which, h{s), is set by the conservation of the total growth probability, 
J dsp{s) = 1. Note that in Eq. (3), two factors of p{s') appear-one corresponds 
to the original probability that a particle lands at p{s'), the other to the potential 
trajectories arriving at s that are blocked by such a particle. 

For a <^ \s — s'\ <^ an, conformal transformation shows that the function 
H{s, s') is given in two dimensions by the simple form^ 

ms.s')= ^ i'WpW (4) 

J^ ds"p{s") 

where the integral in the denominator is the total growth probability between the 
points s and s'. It is convenient to parameterize the interface by this quantity, 
the "growth probability" distance between points z{s), defined by z{s') — z{s) = 
J ds" p{s"). Then our fundamental equation becomes 



dn 



^ h{z)6{z-z'] 



p\z'), (5) 



{z-z'Y 

where a serves as an ultra-violet cutoff to prevent divergence of the integral, and 
h{z) is related to /i(s) and to the function -2(s); its precise form is of no interest to 
us. 

I wish to use this equation to determine the function dx/dlnnb = g{x,y). 
Repeated application of the chain rule yields 



d\nn,~ pIY ' dn dn j ' ^^' 

Consider a branch with probabihty p' and a number of particles n' . We suppose 
that this branch extends from z = to z = p' . Eqs. (5) and (6) imply that if we 
can write p'^{z) on this branch (and by extension, all other branches) as 



P'W = ^/(^/pO, (7) 

n 



where f{z) is a universal function that depends neither upon p' nor upon n', then 
we will be able to write dx/dlnub = g{x, y), with the right-hand side a function of x 
and y alone. Equation (7) is motivated by the fact that p'^{z) must be proportional 
to (p')"^; the dependence on n' is specifically chosen to lead to an n'-independent 
g{x, y). Only if we can find a method of computing an n'-independent f{z) will this 
ansatz be justified. 

Thus the crux of the problem is this "branch envelope" function f{z), which 
represents, with the appropriate normalization, the distribution of growth probabil- 
ity in different regions of a branch. Now in our picture, each branch can be divided 
into two distinct sub-branches, which compete according to the dynamics estab- 
lished by g{x,y). Our central mean-field assumption is that we can compute f{z) 
by averaging the envelope functions f{z) of these sub-branches over the stochastic 
parameter e appropriate to the competition of these two sub-branches. In this way 
we obtain the following equation: 



m-i -M.^S/(^).^i^S/(,^l^ (8) 
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where x{eni,) and y{enb) give the values of x and y along the unstable manifold as 
functions of n^ and the stochastic parameter e. For convenience, we are defining p(e) 
for negative values of e as p(— e) = /o(e), with x{—r]) = 1 — x{r]), y{—rf) = 1 — yivj). 
This leads to the relatively compact expression of Eq. (8). For large n^, this equation 
has a solution independent of n^, which is determined by 



dr]\r] 



u-l 



x'^{r]) 



f 



I ii-x{v)? f 



y{v) \x{v)J i^-yiv)) \^-x{v) 



\- z 



/(^) = 0. (9) 



Since the integrand goes to zero as ?7 -^ oo, we are justified in taking the small e 
form for /9(e). 

Of course, in order to perform this integral, we must have the form of the 
unstable manifold, and thus we must already know g{x, y). We can determine g{x, y) 
from f{z) by simply integrating Eq. (5) over the appropriate intervals. We do not 
integrate over regions exterior to the two competing branches, but only investigate 
the infiuence of the two branches on one another. Skipping some tedious algebra, 
we may express the result as follows. Defining a function t(j{u) by 



t(j{u) 



dz { - 
V^ 



1 



z + u 



/w 



(10) 



we can write 



g{x,y) =x{l -x) 

1 — X 



2-^(oo)-^^^ 

y {^-y)x 



ij 



X 



'(/'(oo) — 



X 



y{l-x^ 



■^ 



1 

1 — X 
X 



X 



(11) 



The reader should note that we have a circular procedure, because g{x, y) is 
determined as a function of f{z) by Eqs. (10) and (11), while f{z) is determined as 
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a function of g{x, y), and in particular by the unstable manifold in the x — y plane 
as determined by g{x, y), by Eq. (9). Thus in practice we are looking for a solution 
of Eq. (9) where the functions x{ri) and y{r]) are implicitly determined by f{z). 

I have numerically obtained the unique solution to Eq. (9) under these condi- 
tions, which is displayed in the inset to Figure l.^'^ This validates our assumption 
regarding the scaling with n' in Eq. (7). The function g{x,y) determined from this 
function has all of the necessary qualitative features; in particular, the fixed point 
at {x,y) = (1/2, 1/2) is unstable and hyperbolic, and the unstable manifold leads 
from this point to stable fixed points at {x,y) = (0,0) and (1, 1), as illustrated in 
Figure 1. Figure 1 also shows numerical results for branch competition. The value 
of the unstable eigenvalue z/ is z/ ^ .6020, implying that D = 1/u ^ 1.661, which is 
within 3% of the standard numerical result D ^ 1.71. 

In addition, this theory automatically agrees with the electrostatic scaling law, 
which states that 

{ ds p{sf) (X n-^ , (12) 



where the integral is over the entire cluster surface. This is equivalent to the more 
usual statement that r(3) = D. In ref. 4, we demonstrated that the multifractal 
exponents a{q) defined by J ds p{s)'^ oc n~^^'^^ can be obtained from the integral 
condition^ ^ 



<"»■»-' 4^^ + ,r ::":„. -1-0- ds) 



°° ^ x{ri)i {l-x{ri))i 

""'"' {yW^^ (1 - 2/(77))- (<?) 

By integrating Eq. (9) from z = to z = 1, one obtains precisely this criterion, with 
q = 3 and a{q) = 1, in agreement with the electrostatic scaling law. Though the 
electrostatic scaling law thus appears in a natural way in this theory, one should 



10 

not say that it is predicted by this theory unless the solution obtained to Eq. (9) 
is stable. It may be that it is necessary to impose the electrostatic scaling law as a 
constraint to insure this stability. ^° 

From Figure 1, it is clear that although in some sense the unstable mani- 
fold that we have calculated is an acceptable average trajectory, the numerically 
obtained trajectories do exhibit some dispersion about this average. This has sig- 
nificant results. The Makarov scaling law predicts that da{q)/dq\q=i = I/Dq,^'^ 
where Dq is the surface fractal dimension (which according to some studies is sig- 
nificantly less than the radius-of-gyration dimension D.)^'^ My result, from Eq. (13), 
is da{q)/dq\q^i ~ 0.71, which is significantly different from the Makarov result. In 
practice, this quantity is quite sensitive to the way in which the unstable manifold 
approaches the stable fixed points at (x, y) = (0, 0) and (1, 1); since the numerical 
trajectories are quite dispersed in this region, I do not expect a good result for the 
Makarov scaling from a one-trajectory theory. However, the theory outlined in this 
letter can be easily generalized to account for the possibilty of trajectory dispersion, 
which may lead to better agreement with the Makarov result. 
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Figure Caption 

1. Trajectories of branch competion in ttie x — y plane. Tlie ligtit solid trajectories 
are numerical results from ref. 4 for specific branch pairs in growing DLA 
clusters. The heavy solid line represents the unstable manifold predicted by 
this letter, which is quite close to the "average" numerical trajectory. The inset 
shows the computed branch envelope function f{z). 



